function diff = chi_fun( param)
    % vacancy filling rate
    diff = (param.lambda./param.lambda_mom).^( - param.epsilon./(1-param.epsilon)); 
    
end

